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}Q . Abstract 

T^lj- I We study in detail the dynamics of a scalar field in thermal bath with symmetry 

^^ ' breaking potential. In particular, we focus on the process of trapping of a scalar field 

at an enhanced symmetry point through the thermal/non-thermal particle produc- 
tion, taking into account the interactions of produced particles with the standard 
model particles. As an explicit example, we revisit the saxion dynamics with an 
K> ' initial amplitude much larger than the Peccei-Quinn scale and show that the saxion 

C^ , trapping phenomenon happens for the most cases and it often leads to thermal infla- 

tion. We also study the saxion dynamics after thermal inflation, and it is shown that 
thermal dissipation effect on the saxion can relax the axion overproduction problein 
from the saxion decay. 
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1 Introduction 

After the discovery of the Higgs boson at the LHC [1], the concept of spontaneous sym- 
metry breaking induced by the vacuum expectation value (VEV) of a scalar field becomes 
more important. Actually, apart from the electroweak symmetry breaking, the Universe 
may have experienced several phase transitions such as the Peccei-Quinn (PQ) phase tran- 
sition p], grand unified theory (GUT) phase transition and so on. 

Most of these phase transitions are thought to be associated with VEVs corresponding 
scalar fields, cj). While the concept is simple, the dynamics of the scalar fields leading to the 
symmetry breaking might be rather complicated and may have significant cosmological 



implications. The usual argument is as follows [3]. At high temperature, sits at the 
origin = due to thermal effective potential since it is an enhanced symmetry point. As 
the temperature decreases, thermal effects become less significant, and finally relaxes to 
the minimum. However, the assumption that starts at the point = may not always 
be justified, because can have very large initial value during and just after inflation. 
Then we need to study dynamics of a scalar in order to see whether or not the symmetry 
is restored (i.e., = 0) in the early Universe. 

If the initial amplitude of the scalar field is displaced from the minimum of the potential, 
the subsequent dynamics of the scalar field is non-trivial. With such an initial condition, 
the scalar field oscillates and dissipates its energy via the interaction with background. If 
the effect of the cosmic expansion is negligible, the scalar field eventually relaxes to an 
equilibrium state, which is determined by the initial energy (and other conserved quanti- 
ties). In such a case, the scalar field is expected to be trapped at the enhanced symmetry 
point if the temperature of the equilibrium state is higher than the critical temperature 
for the phase transition. On the contrary, if the cosmic expansion alone efficiently reduces 
the energy of oscillating scalar field, the scalar field cannot experience the restoration of 
the symmetry and results in the symmetry breaking vacuum. The actual situation lies 
somewhere between these two extreme cases. Therefore, to see whether the scalar field 
experiences the trapping or not, we have to study in detail the dynamics of scalar field and 
compare the time scale of dynamics with that of cosmic expansion. On this respect, the 
phenomenon called non-thermal phase transition [^ and moduli trapping [5] was studied 
previously. 

The complete analysis of dynamics of scalar field in a concrete setup including inter- 
actions with standard model (SM) sector in thermal environment has been lacking. In a 
realistic setup, such a scalar field often interacts with gauge-charged matters and the back- 
ground is filled with thermal plasma after inflation. The primary purpose of this paper is 
to analyze the scalar dynamics with symmetry breaking potential by taking into account 
all possible relevant effects: thermal correction to the scalar potential, dissipative effects 
on the scalar field, non-perturbative particle production and the resulting modification on 
the scalar potential. 

To make our discussion concrete, we will pay particular attention to the dynamics of 
saxion in a class of supersymmetric (SUSY) axion model. We will show that, even if 
the initial amplitude of the saxion is much larger than the PQ scale, the zero-mode of 
the saxion field is eventually trapped at the symmetry enhanced point in large parameter 
space; this is due to the modification of the saxion potential and the dissipation of the 
oscillation. Thus, the thermal inflation [SJ [TJ E] is quite often caused by the saxion. 
We also study the scalar dynamics after phase transition. We find the situation where 
thermal dissipation into SM plasma efficiently reduces the saxion energy density without 
axion overproduction. 

This paper is organized as follows. In Sec. [21 we introduce a basic setup and review 
thermal effects and particle production processes as preparations for the analyses in the 
following sections. In Sec. [31 we study the saxion dynamics in detail in the cases of both 



small and large mixing between PQ and SM quarks. We will study in detail whether 
the saxion is trapped at the origin or not by taking account for thermal and non-thermal 
effects in the following sections (Sec. I3.2l and Sec. I3.3p . In Sec. HI we will discuss the saxion 
dynamics after thermal inflation and its cosmological consequence. Sec. |5] is devoted to 
the conclusion and discussion. 

2 Preliminaries 

2.1 Overview 

Here, we consider the evolution of a complex scalar field which is responsible for the 
spontaneous breaking of a global U{1) symmetry; the (zero-temperature) potential of 
has a minimum at \(j)\ = v y^ 0. Above this scale, the potential of (p, which is lifted by the 
SUSY breaking effect, is assumed to be quadratic: 

r ~mVP for 101 >t;. (2.1) 

In addition, the scalar is assumed to have interactions with scalar Q and fermion Q as 

£ = A^|0ngp + (A;0Qg + h.c.), (2.2) 

with coupling constants A^ and A/. Q and Q are assumed to have SM gauge interactions. 
Having a SUSY model in mind, we call Q and Q as quark and squark, respectively, 
hereafter. 

How does the phase transition proceed? As described in the Introduction, the usual 
arguments are as follows. At high temperature, sits at the origin = due to thermal 
effective potential. In our case, it may be caused by the interaction with Q and Q. As the 
temperature decreases, thermal effects become less significant, and finally relaxes to the 
minimum |0| = v. 

However, it is non-trivial whether actually starts at the point = 0. Let us assume 
that initially, just after inflation, is placed far from the origin: |0| 3> v. This assumption 
is quite natural if m ^ -f^inf, where ifinf denotes the Hubble scale during inflation. An 
immediate consequence is that since Q and Q get large masses from the VEV of 0, they 
may decouple from thermal bath even after the reheating, depending on the reheating 
temperature. Thus we need to investigate the dynamics of in order to see whether 
will be eventually trapped at = or not. 

• If 101 is large enough, Q and Q decouple from thermal bath. Still, however, feels 
thermal effects through the so-called thermal logarithmic potential [9], which may 
tend to stabilize toward the origin. 

• The dissipation is essential to effectively reduce the coherent oscillation energy den- 
sity [TOl O [m HH Uni [m ^5\. in particular, we must take into account thermal 
dissipation effects on the scalar in order to see truly relaxes to the origin. 



• In the limit of low reheating temperature, such thermal effects would be inefficient. 
Even in such a case, because of the large amplitude of 0, non-perturbative particle 
production events happen at each oscillation of (j). Then, the finite number density 
of Q and Q modifies the effective potential of which also tends to stabilize (j) at 
the origin. 

• The produced Q and Q decay or annihilate depending on their interactions to the SM 
particles. If Q has the same gauge quantum numbers as of right-handed down-type 
quarks, for example, we may obtain the following mixing term 

Cra\^ = nQqiH + ]i.c., (2.3) 

where qi and H are SM left-handed quark and Higgs boson, respectively. With 
the above interaction, Q soon decays and generates thermal plasma. Even without 
such mixings, annihilation processes into SM particle also reduce Q and Q, which 
would also generate thermal plasma. Through these processes, coherent oscillation 
effectively lose its energy by producing thermal plasma, which may significantly affect 
the scalar dynamics. 

Thus the actual dynamics of the scalar would be complicated due to the combinations 
of these effects. We will analyze the scalar dynamics in Sec. |3]with an explicit example of 
the saxion in SUSY axion model, but we emphasize that the analysis can be applied to 
general class of scalar fields with symmetry breaking potential. In the rest of this section, 
we briefly summarize the thermal effects and particle production processes. 

2.2 Thermal and non-perturbative effects on scalar dynamics 

When the primordial inflation ends, the scalar field may be displaced far from the true 
vacuum |0| = v. Eventually, it begins to relax towards the equilibrium state which is not 
necessarily identical to the true vacuum because there is a background thermal plasma. 
Even with this initial condition, which is far from the origin of potential, the scalar field 
may be trapped at the origin due to thermal effects since the (s) quarks become massless 
at the origin of scalar potential. In addition, even if thermal effects are not efficient, 
the scalar field can be trapped by the explosive production of squarks via the parametric 
resonance [161 IE] • 

Now let us summarize the effect of background thermal plasma at the beginning of 
scalar coherent oscillation and also the conditions for non-perturbative particle production 
in the presence of thermal plasma [151 IIH] • We will discuss whether the saxion is trapped 
at the origin or not by taking account for thermal and non-thermal effects in the following 
sections (Sec. 13.21 and Sec. 13. 3p . 



2.2.1 Thermal effects at the beginning of oscillation 

Since the (s)quarks are in the thermal plasma if A|0| < T, there might sink at the origin 
of free energy. At the one-loop level, the free energy is given by [19] 

rp4. /»00 



where i represents the species (normalized by one complex scalar or one chiral fermion), 
|i| denotes the statistical property: \i\ = 0, 1 for boson and fermion respectively, and M, 
is the mass of i that depends on |0|. If the mass is smaller than the temperature M, -C T, 
then one can find that this term encodes the so-called ^^ thermal mass" term: 

^i-ioop3 J2 ^T'M^i\<f>\)+ E ^T'MU\<j>\). (2.5) 

igBoson iGFermion 

On the other hand, if M, ^ T, this one-loop contribution rapidly vanishes. Instead, 
higher loop effects dominate the free energy, which arise from the threshold correction to 
the gauge coupling constant g at T. That is, using the fact that the free energy of hot 
plasma has a contribution proportional to g^{T)T^, and that the gauge coupling constant 
at the scale below Mi has a logarithmic dependence on the scale Mj, one finds the free 
energy of 0, so-called ^Hhermal log" potential as [9j: 

Fth-iog = ai^aiT^THn [AVIV^^] (2-6) 

where cl is an order one constant. 

Thus, the free energy that affects the dynamics of scalar condensation can be approx- 
imately parametrized as 



Kh = < (2.7) 
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where cm/l are order one constants. The effective potential for the scalar condensation is 
given by the sum: Ves = V + Vth- 

The scalar condensation begins to oscillate when the Hubble parameter becomes com- 
parable to its effective mass: 

Has ^ max [m, XT (for A0i < T), aT^/^i (for T < A0i)] (2.8) 

with the initial amplitude (pi being assumed to be larger than v. Here we omit the order 
one constants cm/l for simplicity. There are three cases depending on which term domi- 
nates the effective potential: thermal log, thermal mass and zero temperature mass. The 
temperature Tqs at the beginning of oscillation for each cases is summarized as follows [15] : 



The scalar begins to oscillate with thermal log potential if 



< aTu 



M, 



pi 



2/3 



m 



and 0, > (Tr/A)^^^ (aMpi) 



,1/3 



The temperature at the beginning of oscillation is given by 



T 



laM, 



pi 



Tr. 



(2.9) 



(2.10) 



In addition to above two inequalities, the condition Tqs > Tr should be met since 
otherwise the scalar oscillates with the zero-temperature mass term [case (iii)]. 

The scalar (j) begins to oscillate with thermal mass if 

nl/4 



A> 



A> 



m 



TlM^i 



and 4>i < (Tr/A)^/^ (ttMpi)^/^ for AMpi > T, 



R 



(2.11) 



m 



M, 



pi 



for AMpi<Tr. 



The temperature at the onset of oscillation is given by 

[AT^Mpi]'/' for AMpi > Tr 
AMpi for AMpi < Tr. 



T 



(2.12) 



Otherwise, the scalar begins to oscillate with zero temperature mass, and the 
temperature is given by 



-.21 1/4 



T 



[mM^xTl]' for mMpi > T^ 



'mM, 



pi 



for mMpi < T|. 



(2.13) 



2.2.2 Non-perturbative particle production 

The dispersion relations of (s)quarks depend on the scalar field value: 



^qW = y/A2|0(t)P + m,\,Q + k2 



(2.14) 



with the soft (s) quark mass being neglected. Here mefr,Q encodes the finite density cor- 
rection to the real part of dispersion relation: e.g. the contribution from particles in the 
thermal plasma is imprinted as the thermal mass term. In our case, the initial amplitude 
(j)i is much larger than the scalar mass m: X(f)i ^ m, and hence perturbative decay is 
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kinematically forbidden during the most region of oscillation period. Then, the following 
non-perturbative particle production becomes important. 

The non-perturbative particle production occurs when the adiabaticity of (s) quarks is 
broken down: IcJq/wqI ^ 1 [TBI |17]. The condition for non-perturbative production is 
given by 



A0 ^ max 



m, 

m 



(2.15) 



where (f) is the amplitude of the oscillation of (f). Note that even if the back-reaction of 
particle production dominates the dynamics of scalar, the condition [Eq. fl2.15p ] is useful by 
simply replacing m with the effective mass rhes that encodes the finite density correction. 
If Eq. (12.1 5p is met, after the first crossing of the region where the adiabaticity is broken: 
101 < 0NP with 0NP = (m0/A)^/^, the number density of (s)quark suddenly acquires the 
value 

UQ^^,; K^{\m~^f/^ (2.16) 

for one complex scalar or one chiral fermion [17]. Before the first crossing of the region 
where the adiabaticity is broken, we have uq = 0. Thus, from Eq. (l2.15p . we obtain the 
condition for the non-perturbative production: 

A;2 > max [m2,/T2] . (2.17) 



In parameters of our interest |^ the non-perturbative production may occur when the 
scalar begins to oscillate with the zero-temperature mass, and the reheating takes place 
after the beginning of oscillation. Hence, at the first crossing, the latter condition of 
Eq. fl2Tri) implies 

/ 777, 

rR<(A//)J— 0.. (2.18) 

Then there are three cases. 

(A) begins to oscillate with thermal log potential. 

(B) begins to oscillate with zero-temperature potential without non-perturbative par- 
ticle production. 

(C) begins to oscillate with zero-temperature potential with efficient non-perturbative 
particle production. 

We discuss the scalar trapping in these three cases separately. 

In Fig. m we show the contours of constant Tqs on 0j vs. Tr plane; here, we used 
171 = 1 GeV, A = 0.05, and a of the strong interaction. The shaded region (A) corresponds 



"^^ Here we assumed X < g"^. See also footnote |#2 
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Figure 1: The contour plot of Tqs as a function of (pi and Tr,. Here, we have taken m = 1 GeV, 
A = 0.05, and used a of the strong mteraction. In the region (A), the scalar begins to oscillate with 
thermal log potential. In the region (B), the scalar begins to oscillate with zero-temperature mass 
without non-perturbative particle production. In the region (C), the scalar begins to oscillate 
with zero-temperature mass with efficient non-perturbative particle production. 



to the parameters where the scalar begins to oscillate with the thermal log term. On 
the other hand, in the regions (B) and (C), the scalar begins to oscillate with the zero 
temperature mass term. The boundary between (A) and (B) [(B) and (C)] is given by 
Eg. (12. 9 p [Eq. (I2.18P ]. Note that the region where the scalar starts to oscillate with the 
thermal mass term lies outside of this figure: higher Tr and smaller 0j. We will not 
consider this case in the following. 

The subsequent evolution of scalar /plasma system after the first crossing crucially 
depends on the setup. There are two cases: 

(i). If the produced particles cannot decay faster than the oscillation period, the para- 
metric resonance can occur due to the induced emission effect of previously produced 
particles. Hence, the number density grows exponentially and then the back-reaction 
soon dominates the effective mass of (f). In addition, because of the largeness of the 
number density of Q in the background, the linear potential of cj) is induced. [See 
Eq. (I2.23p .] This case will be extensively studied in Sec. 13.21 

(ii). On the other hand, if the non-perturbatively produced (s)quarks can decay into other 
light particles much faster than the oscillation period of (^, then the phenomena like 
instant preheating pO] take place. In this case, the linear potential generated by 



the produced particles is insignificant since the produced particles soon decay. See 
Sec. 13.31 for further discussion. 

First, let us consider the case (i). After the passage of non-adiabatic region \(p\ < 0np, 
the produced particles become heavy due to large field value |0(t)| and their decay rate 
become large correspondingly. If the decay rate Fj,^ ~ K'^\\(f){t)\ [with k being the mixing 
parameter given in Eq. ( 12. 3p ] reaches T^e^{tdec)tdec ~ 1 with tdec ^ m~^, the produced 
particles completely decay into SM radiation well before the scalar moves back to the 
origin. The condition for the case (i), which requires the stability of Q for the time scale 
of our interest, is given by 



/ m 

V ^0 



(2.19) 



If the above condition is satisfied, the number density of (s)quarks grows explosively due 
to the parametric resonance as (p crosses the region where the adiabaticity is broken. In 
our set up, the condition Eq. fl2.15p implies 



k^ > max 



2 ^2rp2 Jl'^Q \2^Q 

off , 9 -t ,9 -: — , A - — 



^e^^9 



fC* 



(2.20) 



where rries represents the effective mass for that includes the back-reaction of particle 
creation and it is estimated as 



r/igg ~ max 



2 x^Q 



(2.21) 



Eventually, the parametric resonance stage finishes when the inequality Eq. (I2.15P is sat- 
urated. This condition can be rewritten as 



rriefi <^ min 



k 



X' 



h 

*5 2 * 



fc* 



(2.22) 



In the following, we consider the case: A < (^^ for simplicityl^^land hence c = g^ jy? ■ In 
addition, in the case (i), the linear potential of 0, which is obtained as [17] 



Vx 



linear 



A^I0P(IQH~A^ 



d^k 



(27r)3u;Q 



JQik) ~ A|0|nQ for 101 > 



NP, 



#3 



(2.23) 



(Here, 



is generated because sizable amount of Q exists irrespective of the field value of i 

fq denotes the distribution function of Q.) Such a linear potential plays very important 

role in studying the evolution of 0, as we will see in the next section. 



'^^ For large A, (/> may be more likely to be trapped at the origin since the region (C) where the non- 
perturbative production occurs becomes larger [See Eq. ()2.18p ]. 

"^^^Here we used the fact that the k integration is dominated at A; ~ fc* and ujq ^ A|0|. For the case of 
|0| < (/)NP, WQ ~ fc*, and hence the potential becomes quadratic in |(/>p. 



Next, let us consider the case (ii), which corresponds to the case where the condition 
f l2.19p is not satisfied. In this case, the produced particles soon decay at |0(tdec)| ~ 
[m(j)/{K^\)Y''^ <^ (p. Taking an oscillation time average, one finds the effective dissipation 
rate in this case [15j: 

rf^^^iVdof x^ (2.24) 

where iVd.o.f. stands for the degree of freedom normalized by one complex scalar or one 
chiral fermion. In addition, in case (ii), the parametric resonance does not occur and the 
linear potential generated by the produced particles is insignificant. 

3 Saxion Dynamics 

The strong CP problem is one of the serious theoretical problems in the SM. The PQ 
mechanism is a promising solution to the strong CP problem [2]. It utilizes a global U(l) 
symmetry, called PQ symmetry, which is spontaneously broken at the scale fa- The axion 
appears as a consequence of spontaneous breakdown of the PQ symmetry and it plays an 
important role in cosmology and phenomenology [2T| 122] . 

In SUSY extensions of the axion model, the early Universe cosmology becomes much 
more non-trivial because of the presence of fiat direction in the scalar potential, called 
saxion [231 EH EHl [26]. The cosmological evolution of the saxion field depends on the 
detailed structure of the saxion potential. Among various possibilities (see Ref. [22] for a 
partial list of the stabilization mechanism), in this paper we focus on the case where there 
is only one (complex) PQ scalar field having VEV [271 |28l [291 130] . 

The potential of the saxion field is lifted only by the effect of SUSY breaking and 
is very flat for |0| ^ /„. Then, without a tuning, the saxion is expected to be initially 
placed at the origin (|0| = 0) or around the Planck scale (|0| ~ Mpi), depending on 
whether the saxion receives positive or negative Hubble mass correction during inflation. 
In the former case, the saxion likely causes thermal inflation [6l [71 [8] as extensively studied 
in Refs. [311 [321 [331 EH [35] . Once thermal inflation happens, the saxion must successfully 
reheat the Universe before the big-bang nucleosynthesis (BBN) begins. Unfortunately, the 
saxion often dominantly decays into the axion pair, which would invalidate the thermal 
inflation scenario. (However, it is possible that the saxion dominantly decays into the 
Higgs sector in the DFSZ model |36j.) 

On the other hand, the case of large initial amplitude of the saxion, \(p\ ^ /„, has not 
been studied in detail so far. We mainly focus on this case and investigate its cosmological 
effects. The dynamics of the saxion in such a case is rather complicated due to two effects: 
thermal effects and non-perturbative particle production]^ As will become clear, both 
effects tend to temporary stabilize the saxion around the origin, |0| = 0. We will show 



"^"^ Thermal effects on the saxion dynamics were considered in Refs. [371 [H] in a different class of PQ 
model. 
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Table 1: Charge assignments 
^ Q Q 

U(l)pQ +1 -1 

SU(3)c 1 3 3 



that, even if we start with a large saxion initial amplitude, it will be eventually trapped 
at the origin due to these effects and thermal inflation takes place for the most reasonable 
parameter choices. Cosmological issues related to this will be discussed in Sec. HI 

3.1 Setup 

Amongst several ways to stabilize the saxion, we focus on the models with only one PQ 
scalar whose potential is stabilized by the quadratic term from SUSY breaking [271 l2Hl ISH 

m\. 

To capture the essential features of saxion dynamics within these models, we approxi- 
mate the saxion potential as 

-m2(|0|2-M2) for 101 <M 

V ~ { -m^M^ log^(|0|/M) for M < |0| < /„ (3.1) 

-m'm'-f!) for /,<|0|. 

where fa is the PQ scale, rus and m are mass parameters terms from SUSY breaking effect. 
(Note that we have the relation mfa — rrisM.) In addition, M is the messenger scale in 
the framework of gauge-mediation. (For the cases of gravity- and anomaly-mediation, see 
the discussion below.) The relation between rris and m depends on the mechanism of 
SUSY breaking. In the models based on gauge-mediation |27l IMl EHl SO] , "^ is expected 
to be of order the gravitino mass. Furthermore, the negative soft mass of saxion potential 
at its origin comes from the three loop effect via the Yukawa coupling A, and hence it 
is related with the squark mass as wigoft oc nis/X (with rrisoit being the sfermion masses 
in SUSY SM.) As a result, there can be a hierarchy between M and /„. In the gravity- 
or anomaly-mediation models, the renormalization group effect may stabilize the saxion 
potential ^ i29j. If so, m ^ rus (and M c^i fa). 

We consider the SUSY extension of the hadronic axion model [41]. We introduce the 
following superpotential: 

W = X^QQ, (3.2) 
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where Q and Q are identified as the PQ quark superfieldqi^ charged under both SU(3)c and 
U(l)pQ so that it has U(1)pq-SU(3)c-SU(3)c anomaly and solves the strong CP problem. 

The charge assignments are summarized in Table [Tl^i (Here, the PQ charges of Q and Q 
are arranged so that the mixings with the SM fermions are possible.) 

With the PQ charges given in Table [H there may exist the following terms in the 
superpotential 

W^mix = t^QqiBd, (3.3) 

which induce the mixing with the SM particles. (Here, all the particles in the SUSY SM 
are vanishing PQ charge.) If these mixing terms exist, PQ quarks decay into SM fields. As 
will become clear, the dynamics of saxion significantly depends on whether PQ (s)quarks 
are stable or not. In this paper we consider two cases separately: small mixing limit 
(Sec. 13. 2p and large mixing limit (Sec. 13. 3p . 



3.2 Saxion dynamics: small mixing case 

As mentioned at the beginning of Sec. 12.21 the saxion may be trapped at the origin due to 
the thermal/non-thermal effects even though its initial field value is displaced far from the 
origin. Let us discuss thermal and non-thermal effects in turn and derive the condition 
for saxion trapping. In this subsection we focus on the case of small mixing limit between 
PQ and SM quarks, and hence PQ (s) quarks are considered as stable particles. 

3.2.1 Case (A): Thermal trapping 

First, we consider the case (A). If the negative mass of saxion potential at the origin 
vanishes due to the thermal effects when the Hubble parameter becomes comparable to the 
dissipation rate of saxion into radiation, the saxion is trapped at the origin. In our setup, 
the dissipation rate is expected to become larger than the Hubble parameter soon after the 
onset of oscillation because the temperature Tqs is much larger than m and the coupling A 
is not so small. Hence, we derive the condition for saxion trapping with the approximation 
that the saxion immediately dissipates its energy after its onset of oscillation. Note that 
there is no explosive particle production due to the non-perturbative effect for A < a. 

Let us consider the first crossing of = 0. The time scale St, during which the saxion 
passes through |0| < 0^ = T^s/X, is estimated as 6t ~ l/(AaTos). The abundance of PQ 
(s)quarks that is produced during this time interval is given by nq ~ n^Y^{o'prodi^)St, where 
((^prod"^) ~ 9*/Tos denotes the PQ (s)quark pair production cross section and rith ~ T^^ is 



"^^ Symbols with hats represent the superfield. The corresponding scalar component is represented by 
symbols without hats. 

"^^There can also be PQ leptons, L and L, which are combined with PQ quarks to form SU(5) fundamen- 
tal and anti-fundamental representations as {Q,L) and {Q,L). Even with those particles, the following 
arguments are unchanged. 
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the number density of particles in thermal bath. Thus the number density after the first 
passage for |0| < 0^ is evaluated as 

nh ^ Ti, (3.4) 

as long as A < a. (Here, the superscript "s" is for the saturated value.) After the 
passage, the saxion may climb up the linear potential Ves — ^T^s\^\ ^^^ \^\ > 4>c, and 
eventually reach its maximum value in 5tmax ~ a /{XT). After the decoupling of PQ 
squarks at |0| > (pc, the annihilation of PQ squarks may still be efficient. In this case, 
their abundance reduces to uq ~ {{crannv)Stjas,x)^^ ■ The result for |0| > 0c is 



Uq ~ min 



1.4 



X r„t ^ rf,T„t. (3.5) 



Recalling that we are considering the case where PQ (s)quarks are (quasi-) stable, we 



obtain the effective potential of the saxionj^^ 



— ^|0P for 101 < 0, 
/eff '^ s (Pc (3.6) 

An' 101 for 101 > 0c. 
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Then, the conservative condition for saxion trapping is summarized as follows: 

• For Tos > AM, the saxion can be stabilized at the origin solely by the effective mass 
term. The condition is nothing but rri^ < XnQ/cf)^- This condition is rewritten as 



1/2 /l\i/' fmlX^ ^ -*- ^'/' 



For Tos < AM, the combination of the effective mass and linear term stabilizes the 
saxion at the origin. The condition is given by rn^M'^ < Xn^qM . This is rewritten as 

\ Aa\ J \aivip\j 

/ 1 \ /™ /\ \ 1/3 / ^f \ 1/3 / X \ 1/2 

. (3.8) 

Note that these conditions are met for most parameters of our interest if the saxion oscil- 
lates with the thermal effects (recall that Eq.( ]2.9p : aT-^ > (pi^Jm/Mpx in the case (A)). 
The saxion oscillation around the effective potential soon decays because the effective 
dissipation rate, r\ , is large: P^ ~ X^aT. Eventually, the saxion is trapped at the 
origin. 




^^T4ff can be evaluated as Viinoar given in Eq. (|2.23p . Note that, in the present case, the k integration 
is dominated at fc ~ T. 

"i^^Throughout this paper, we impose the foUowing conservative condition; we require that, for the 
trapping at the origin, the origin be the absolute minimum of the thermal potential. We do not consider 
the case where there are other minima than the origin. As we will see, even with this conservative criteria, 
the saxion is likely trapped for most parameters of our interest. 
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3.2.2 Case (B): Quasi-thermal trapping 

Next, we consider the case (B). In this case, the thermal effects on the saxion potential 
is not important at the beginning of oscillation. Still, however, the situation is similar 
to the case (A) studied above. When the saxion passes through the origin, the thermal 
bath creates PQ (s)quarks. To see this, let us estimate the time scale 5t, during which the 
saxion passes through the region |0| < 0^ = T^s/X- It is estimated as 6t ~ Tos/k^ where 
/c* = yJXuKpi. The number density of PQ (s)quarks produced during this time interval is 
given by hq ~ nl^{(y^^o^v)5t, and it is estimated as 



Un — nim 



1,- 
e 



X Tl ^ d^Ti, (3.9) 



where e = kl/{g'^T^J and 4 = {l + e/g'^)~^. Note that A < e < 1 in the case (B). If e > g'^, 
PQ (s)quarks does not reach thermal equilibrium during 6t and hence the annihilation 
of PQ (s)quarks is neglected after the passage of |0| ~ 0c- If e < g"^, PQ (s)quarks are 
thermally populated and their pair annihilation takes place after the saxion climbs up the 
potential around |0| ~ 0c- As a result, the number density of PQ (s)quark is given by 

n^Q ^ d'Ji, (3.10) 

where d'^ encodes the model dependent numerical factor: d'^ ~ g'^ /e for e > (7^, and 
d'^ ~ A^/(ea^) for e < g^ and e^a^ < X^ < ea^. The resulting effective potential of the 
saxion is same as (13. 6p . 

Then, the condition for saxion trapping is summarized as follows: 

• For Tos > AM, the saxion can be stabilized at the origin solely by the effective mass 
term. The condition is nothing but m'^ < Xrin/cpc- This condition is rewritten as 

• For Tos < AM, the combination of the effective mass and linear term stabilizes the 
saxion at the origin. The condition is given by m^M"^ < XrigM. This is rewritten as 



Xd'J ^/M~^ 

^ 10-GeV<-/3 (^'^' f^^y^' f^^V' (3-12) 

VlTeVy VlGeV/ VlO^GeVy ^ ^ 

for mMpi > T^. Note that these conditions are met for most parameters of our interest. 
Similar to the case (A), the saxion oscillation around the effective potential soon loses its 
energy due to the dissipation and is trapped at the origin. 
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3.2.3 Case (C): Non-thermal trapping 

Finally we consider the case (C) (fc* > gTos)- Even if the above thermal effects are not 
efficient, the saxion can be trapped at the origin due to the explosive production of PQ 
squarks when the saxion passes through its origin. Let us see this phenomenon, which was 
dubbed as "non-thermal phase transition" 0]. 

As discussed in Sec. I2.2.2[ PQ squarks are efficiently produced for k^ > gT^s for each 
passage of \(f)\ = and its number density exponentially grows until the condition k^ ~ 
TTigg Q is saturated. After the saturation, the amplitude of is given by 



X 



(3.13) 



[For the definition of c, see Eq. f l2.22p .] Using this, the number density of PQ squark is 
given by uq ~ A0g/c^. The saxion effective potential is given by 

Kff - H'^"^^' ^'' l'^' < '^^^ (3.14) 

(^A/iqI^I for 101 > 0NP, 

where 0np = k^/\ ~ (ps/c and m^^ = X^nqjk^ ~ X'^(j)'l/c^. Here k^ = \rhcs(ps with 



K^2 



m 



gfj = XnQ/(f) ~ A^0g/c^. In this case, the squark pair annihilation rate is given by 
Pann.Q = {cF a.nnv) TiQ . Siucc it is dominated by (pit) ~ 0NP, the averaged annihilation rate 
is estimated as 

Eann.Q = ra„n,Q(0Np)^ = YTTI TT*^^ ~ "^^^12^ 2' (^.15) 

Ac* 0NP Xc^ A^ + g'^ 
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Thus the pair annihilation of squarks is not efficient in the time scale of saxion oscillation 
Then, the condition for saxion trapping is summarized as follows: 

• For 0NP > M, the saxion can be stabilized at the origin solely by the mass term. 
The condition is nothing but m^ < rries- This condition is rewritten as 

,.>^.,0»GeVa(^)'(^)-. ,3..) 



"^^ The annihilation of PQ squarks may become efficient shghtly before the condition fc* ^ rrics.Q is 
saturated, since g is not small and the rate of parametric resonance fj-rhcs is smaller than mcS- Even if 
this is the case, the conclusion is almost unchanged: Actually, in this case, the number density of squark 
stops growing slightly before the saturation and the saxion loses its energy in each oscillation via the 
annihilation. The annihilation of PQ squarks produces background plasma. Thus, it is considered that 
the saxion oscillation soon disappears within a few Hubble time and that it is eventually trapped at the 
origin. [See the following discussion in the text below Eq. (I3.2ip .] 
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• For 0NP < M, the combination of mass and linear terms stabilizes the saxion at the 
origin. The condition is given by ml < Xuq/M. This is rewritten as 

'^^^ AV3^ -10 GeVc/A/ ^^^ j (^^^^j l^lO^GivJ " 

(3.17) 

These conditions are likely satisfied for parameters of our interest. 

Here we study the subsequent evolution of the system after the saxion trapping. Af- 
ter the saturation, the saxion loses its energy through the squark pair annihilation into 



radiation j^^i The energy loss in one oscillation of the saxion is estimated as 



„4 4 _ 

^^(0-ann) ^ ruq {a ^nnv) ti^q At ~ -^m^sng ~ —0^, (3.18) 

where At = rh^g using the fact that the energy loss is dominated around ~ 0. We can 
define the effective dissipation rate of the saxion through 

(Q~ann) 



.(Q-ann) _ ^1^^ ^ ^ 9_ 



f W-ann) ^ ^:^_ ^^^ ^ ^^_ ^g^^g^ 



Since this is much larger than the Hubble expansion rate (ifos ~ ^ with ~ 0s, the 
produced PQ squarks efficiently annihilate and correspondingly the saxion amplitude de- 
creases. By solving the equation 



we find 



(3.20) 

-1 -1 
</'W= rs' + -^At-U) . (3.21) 



■p(Q~ann) _ 

+ ^4, P<t>- 


= 0, 


V-'^L^'- 


-U) 



Thus it is seen that the amplitude decays in a time scale ((5t)ann ~ •^c^/(5''^0s), which 
is much shorter than the Hubble time scale. The amplitude follows Eq. (I3.2ip as long 
as the condition for the non-perturbative particle production is saturated (fc* ~ w^efr.Q)- 
The annihilation of PQ squarks produces background plasma, which would give a large 
effective mass to the PQ squarks and terminate the parametric resonance. Actually, after 
a time scale ((5t)ann5 the radiation with a "would be" temperature T^"*-*^-) ~ p/ ~ y/rrKpi 
is produced. (Note that one has T^^-^-) ~ K/X^^^ > {g/X^/^)T > T.) After that, the 
produced plasma soon thermalizes, and the condition for the non-perturbative particle 



'^^^ Since the produced squarks are highly correlated and have non-perturbatively large occupation num- 
ber, the quasi-particle treatment may not be justified soon after the end of preheating, and other processes 
including in-medium-bremsstrahlung and multi-annihilation may also contribute. Here, however, we sim- 
ply estimate the energy transportation of saxion to the radiation by the pair annihilation of quasi-particle 
squarks. 
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production is saturated by the effective mass term for the PQ squark, as kl ~ g^T"^- Similar 
to the above estimates, the saxion also efficiently dissipates its energy into radiation in 
this case. Therefore it is considered that the saxion oscillation soon disappears within 
a few Hubble time around the steep effective potential. Since the temperature is much 
larger than the soft mass scale rris, the saxion is eventually trapped by thermal mass term 
even if the reheating temperature is extremely low. 

After all, for all the cases (A)-(C), the saxion is likely trapped for most parameters of 
our interest. 

3.3 Saxion dynamics: large mixing case 

Now let us study the case of large mixing between PQ and SM quarks, so that PQ (s)quarks 
are unstable. Below we analyze the dynamics of the saxion in the cases (A) - (C) separately 
as in the small mixing case. 

3.3.1 Case (A) 

In the case (A), the saxion begins to oscillate with thermal log potential. As explained 
in Sec. 13.2. ![ the time scale required for the saxion to pass through |0| < 0c = T^s/X is 
estimated as 6t ~ l/(AaTos). The number density PQ (s)quarks produced from thermal 
plasma in this time scale is given by uq ~ nl^{ap^odv)St ~ h?T^^/\. (In the following, 
we assume k ~ 1.) Thus PQ (s)quarks are expected to be thermally populated. After 
the passage, the saxion may climb up the linear potential V^g = AT^gl^l for |0| > 0c, and 
eventually reach its maximum value in (5tmax ~ a /{XT). If the decay time is longer than 
this time scale, the saxion moves back to its origin by the linear potential and consequently 
oscillates with the linear potential. Otherwise, the produced (s) quarks soon decay and the 
saxion oscillates with the thermal log potential. This condition is given by k^ > A/a and 
we concentrate on this case unless otherwise stated. 

Note that the effective dissipation rate of the saxion in the regime of thermal log 
potential is given by T\ ''^ ~ ha^T^ /{cjxp^ and hence we obtain T , /fn^fi ~ Aa. It means 
that the saxion oscillation soon dissipates its energy once it starts to oscillate around the 
origin unless A is very small. This is because the effective dissipation rate decrease more 
slowly than the Hubble expansion rate. As a result, the saxion is trapped at the origin. 

Thus, the condition for the saxion trapping is summarized as follows. 

• For Tos > AM, the saxion can be stabilized at the origin solely by the thermal mass 
term for rrig < ATos- This condition is rewritten as 

-« > T [AT ^ ^°'«^^ iM (^) (lo4v)''- '-^) 

• For Tos < XM, the combination of the thermal mass and the linear potential/the 
thermal log stabilizes the saxion at the origin. The condition is given by rris < 
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aT'^jM. This is rewritten as 

-^ > {mT -^ "'"^ (^) (^'' (^''- '^■"' 

Here we have omitted the case A/a < k^ since it is the same as Eq. (13. 8 p for d\ = 1. 

3.3.2 Case (B) 

In the case (B), the saxion begins to oscillate with zero-temperature mass. As explained 
in Sec. I3.2.2[ the time interval required for the saxion to pass through the region \(j)\ < 
0c = Tos/\ is estimated as 6t ~ Tos/k^. The number density PQ (s)quarks produced from 
thermal plasma in this time period is given by uq ~ n1^{ap^odv)St ~ {k^ aT'^^/ k1)T^^. By 
noting (1/A^ >)a.T'^^/k1 > 1 in the case (B), we find that PQ (s)quarks are expected to be 
thermally populated in this time interval. For most parameters of our interest, the linear 
potential may soon decay, an d he nce the effective potential for the saxion near the origin 
is given by thermal potentialj^lj 

The effective dissipation rate of the saxion in the regime is estimated as V , ~ XaT"^ /(p 

and hence we obtain T\.^^' /m ~ XaT^^/{m(j)). It leads to A^ < P^ ' /m < 1. Thus the 

saxion oscillation soon dissipates its energy unless A is very small since T, "^ decreases 
more slowly than the Hubble rate. 

In summary, the condition for the saxion trapping is as follows: 

• For Tos > AM, the saxion can be stabilized at the origin solely by the effective mass 
term for rris < XT^s- This condition is rewritten as 

_ 1 m^ o /ms/AA / m \-i/2 

• For Tos < AM, the combination of the effective mass and thermal log stabilizes the 
saxion at the origin. The condition is given by rris < aT^^/M. This is rewritten as 

3.3.3 Case (C) 

In the case (C), the non-perturbative production occurs at the crossing of ~ 0. As 
explained in Sec. 12.2.21 the number density of PQ (s) quarks produced in this way is given 



#^^ In fact, the linear potential can decay if 1 < ^^^^St with T^^^ ^ k'^X\(J)\ and dt ^ m(f)/{XT^). Such a 
condition is likely satisfied because T^^^^i ^ {k? /\){m(f>/T'^) > K^a/X. 
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by uq ~ k^/{A7r^). The produced PQ (s)quarks decay at |0(tdec)| ~ [m0/A/t^]^/^. From 
this, we can estimate the effective dissipation rate of the saxion as [15] 



Due to this effect, thermal plasma with temperature T ~ ^yjjrupjjg^ is produced within 
one oscillations of the saxion. This process continues until the condition for the non- 
perturbative production is violated. Actually the c ondi tion is soon violated after a few 
Hubble time after the oscillation, that is, gTjssp ~ k l*^^ [See Eq. fl2.17p ]. Then, thermal 
effects tend to stabilize the saxion at the origin, similar to the case (A) and (B). 

By simply assuming that the temperature of the plasma at the end of non-perturbative 
particle production is given by Tnp ~ y/Xm^i/g, the condition for the saxion trapping is 
summarized as follows. 

• For Tnp > AM, the saxion can be stabilized at the origin solely by the effective mass 
term for nis/X < Tnp ~ y/Xmcpi/g. 

• For Tnp < AM, the combination of the effective mass and thermal log stabilizes the 
saxion at the origin. The condition is given by mg < aT^p/M. This condition is 
rewritten as 47r/a < A0i. 

These conditions are likely satisfied in the most parameters of our interest. Then, as in the 
case of (A) and (B), the saxion soon dissipates its energy unless A is quite small, and sits 
around its origin. Once the saxion is trapped at the origin, it potentially causes thermal 
inflation, as studied in the next section. 

3.3.4 Comment on the case of non-trapped saxion 

As studied above, it is possible that the saxion oscillates with an amplitude much larger 
than the PQ scale without trapping at the origin in the limit of small A and Tr. Then it 
may be non-trivial in which minima the saxion relaxes and whether axionic domain walls 
forms or not. We briefly see what happens in this case. 

Let us suppose that the saxion oscillates in the real axis in with an amplitude and 
define the axion as a = A/2Im0. First note that the mass of a changes from +rn? to —m'^ 
around |0| ^ fa and further it develops to —wig at \(f)\ < M. In one oscillation of the 
saxion, the fluctuation along the axion direction develops due to the tachyonic instability 
as 6a/a ~ /a/0 at |0| < fa- Therefore, the initial axion fluctuation sourced by inflationary 
quantum fluctuations significantly develops at ~ /„ and the motion of is expected to 
become chaotic in the complex plane. This may lead to the formation of axionic domain 
walls, even in the case of no explicit symmetry restoration. 



'^^^ The interaction given in Eq. (13. 3p also affects the thermal mass of Q. Even with such a contribution, 
the present discussion is unchanged because we are considering the case of 5 ~ k ~ 1. 
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A further comment is in order. It was pointed out that the parametric resonant decay 
of the saxion into the axion leads to the nonthermal phase transition |l2l HH HI] . Note that 
these results were based on the potential V oc (|0p — faY- In our setup, the interaction 
term between the saxion and axion exists only in the logarithmic term and the parametric 
resonant enhancement of the axion would be much less efficient. Actually, the condition 
for the broad resonance is only marginally satisfied in our case. In order to make a definite 
conclusion whether the parametric resonance of the axion is efficient or not, we may need 
lattice calculations, which is beyond the scope of our analysis. However, regardless of the 
efficiency of the parametric resonance, we expect the formation of axionic domain walls as 
explained above. 

4 Saxion Dynamics after Phase Transition 

4.1 Condition for thermal inflation 

In the previous section, we saw that the saxion is trapped at the origin = for almost 
all the cases even if the initial position of the saxion is displaced far from the origin. Now 
let us examine whether it causes thermal inflation or not. 

After the trapping, the saxion begins to roll down the potential toward the minimum 
|0| = fa at the temperature T = TpT — rris/X. The condition for thermal inflation by the 
saxion to occur is Prad,inf(^ = ^pt) < m'^M'^^c::^ ""^^/a) where Prad.inf denotes the energy 
density of radiation or the inflaton oscillation, whichever is dominant. 

First, let us consider the low- reheating temperature case: TpT > Tr, so that the 
condition is given by SifprpMpj < m'^M'^. Then, the radiation component in this case may 
contain two contributions: one from the dilute plasma due to the inflaton decay and the 
other from the dissipation of the saxion. The former (latter) is denoted by pr {pf )■ If 
pr dominates the radiation energy density at T = TpT, it is found that thermal inflation 
takes place if 

On the other hand, if pr dominates the radiation energy density at T = Tpx, by noting 
that Pr (Tpt) — m?'(t)1{H-pT:/m)^/^^ the condition is given by 

;. > ^M!!^ . io» GeV f i!^ V'Vi^ V'7^ V^ (4.2) 



Second, let us consider the radiation dominant case: Tpx < Tr. Then the condition is 
given by Tprp < m'^M'^. This condition is rewritten as 

;„>i(I^)^.,«Gev(i|^)(^)\ ,4.3) 
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If either of these conditions are met, thermal inflation is caused by the saxion. Actu- 
ally, unless the reheating temperature is very low or the mass m is small, thermal inflation 
likely takes place. Otherwise, the saxion exits the thermal trapping potential before it 
begins to dominate the Universe and the cosmological problems associated with the axion 
overproduction from the saxion decay, as studied in detail in the next subsection, is alle- 
viated. Even in such a case, however, the saxion may eventually dominate the Universe 
before it decays. In the following, we mainly consider the case of thermal inflation where 
the problem is most prominent. 

4.2 Saxion dynamics after phase transition 

Now let us see the cosmological evolution after the saxion begins to roll down from the 
high-temperature minimum = toward the true minimum |0| = /„. For clarity, we 
assume Tr > TpT hereafter; this is the case in the most of parameters of our interest. We 
divide the saxion dynamics into two stages. At the first stage, the lower edge of the saxion 
oscillation field value is smaller than M. At the second stage, the lower edge of the saxion 
oscillation field value is larger than M. 

4.2.1 First stage 

We consider the dynamics of the saxion after it exits the thermal trapping potential at 
T = TpT ~ rris/X. Let us define the time t so that t = [T(t = 0) = Tpx] as the initial 
condition and reaches |0| = M{fa) at t = tMitf)- Then, we can evaluate the dissipative 
effect on the saxion during its one oscillation around the minimum, using the fact that the 
background temperature is (almost) unchanged for the time scale of one oscillation. 

• t <tM 

The evolution of the saxion is given by |0(t)| = (/)oexp(mst) where 0o is expected 
to be of order T (although its precise value is not important). Thus Im is given by 
tu = "^^^ log(M/0o). The energy loss of saxion due to the dissipation effect when 



it passes through the field value around (j) is given bjl^^ 



,(dis) 



*M 1 . A,^,2t^3 



1 -2 ha^T^ , 2^o , fM\ 



■^ Jo<t<tM Jq 2 (p^ \(poJ 

where fea^ = 9a^ / {128n^ In a'^) ~ 3 x lO^^ ^5]. 

tM <t <tf 

After the passage of \(j)\ = M, the saxion "speed" (/)(~ rrisM) remains nearly constant 



'^-'^^The temperature T ~ mg/A is slightly larger than the mass scale Wg, hence the use of dissipation rate 
in thermal background is marginally justified. 
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until it reaches |0| = fa- Therefore, we have \(p(t)\ = M[H-ms(t— t^)] at this regime. 
We can evaluate the energy loss in a similar way as 



Ap, 



(dis) 



tM<t<tf 



dt -0^ X — ba^T^m^ 

tM 2 02 



(4.5) 



tf<t 

Finally, the saxion climbs up the potential at |0| > fa sXt = tf ^^ fa/irrisM) ~ 1/m. 
The dynamics of saxion at |0| > /a is described by the harmonic oscillation form 
with a frequency m. Parametrically, we have |0| ~ fa and |0| ~ mfa{^ rrisM). Thus 
we obtain 



Ap, 



(dis) 



<></ 



,1-2 fea^T^ , 2^^ 

dt-6^ X ~ ba^T^m. 

2^ 02 



(4.6) 



This is smaller than the energy losses in the regime of t < t /• for m < rris. 



From the above estimates, it is found that the saxion energy loss due to the dissipation 
effect in its one oscillation is given by 



Ap, 



(dis) 



ba^T'^rris log ( -— I • 



We define the effective dissipation rate averaged over the saxion oscillation period: 



(4.7) 



.(dis,l) 



m 



Ap, 



(dis) 



At=m~ 



b'a^T^ ms 
fa m 



(4.8) 



where b' = 61og(M/0o). 

Now let us show the evolution of the radiation and saxion energy density at a few 
Hubble time after the thermal inflation. The radiation energy density evolves as 



(dis,l). 



Prad + 4i/prad = Tl '"' '{T)p^. 



(4.9) 



If r^ "^' ^ -f^PT, the significant fraction of the saxion energy density soon goes into 
radiation and the resulting radiation temperature becomes T = Tti ~ \f^nfa. In the 
opposite limit, after a few Hubble time, the radiation temperature becomes 



where 



/ 



T = max [1, /] TpT = FTp^, 

lO^GeVA 



b'a'X 



M, 



pi 



Ja 



10^ (b'a^X) 



Ja / 



(4.10) 



(4.11) 
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by noting that T ~ const from (14.91) . Therefore, we obtain 

T = min [FTpT, Tti] , (4.12) 

at the end of first stage after thermal inflation. 

By comparing the dissipation rate with the Hubble rate at the end of thermal inflation, 
we obtain 

£ ^ -v.^ . .. . ,....., {^f {l^y (H!|^)' . 

(4.13) 
The ratio (I4.13P can be much larger than unity, meaning that the saxion can dissipate most 
of its energy within one Hubble time after the thermal inflation ends and the radiation 
with a temperature T = Tti is produced. Otherwise, the saxion coherent oscillation soon 
dominates the Universe and it decays into the axion pair, as we will see below, leading to 
unsuccessful reheating. Therefore we must have T\^^' / H-pt ^ 1- 

However, even if r^'^''^'^) is initially much larger than -ffpx, all the energy of the sax- 
ion oscillation is not dissipated away; this is because the dissipation rate is significantly 
suppressed once the smallest value of during one oscillation becomes larger than M. 
Thus, the coherent oscillation around |0| = fa-, whose energy density is the same order of 
the initial total saxion energy density {p^ ~ rn^M'^ ~ rn^f^), still remains after such a 
dissipation. Then we need to consider the second stage of saxion reheating. 

4.2.2 Second stage 

At the second stage, the saxion oscillates around the minimum |0| = /„. Here we assume 
that the condition T\ ^^' ^ Hp^ is satisfied and hence T = Tti. The thermal dissipation 
rate at this epoch is given by 

r?"' = ^. (4.14) 

J a 

By comparing it with the Hubble rate at the end of thermal inflation, we obtain 



(iGiv) [^^) • ^'-''^ 



^p7^^F '°HolMiGivJ [-l^J 

Here we have substituted T = Tpi. 

At the second stage, as well as the dissipation process due to the thermal bath, there 
is another important process which reduces the energy density of the saxion oscillation. 
In the present analysis, we consider the case where the saxion dominantly decays into the 
axion pair. Then, the decay rate of the saxion is estimated as 

1 in 

r.-- = ^-^. (4.16) 

23 



The decay must happen before the saxion again conies to dominate the Universe for the 
successful thermal history. Note that the ratio between Tfp^2a and T\ ^^' is given by 

^X^^ = 6A.ba^ ( ^Y\ (4.17) 

and it is much larger than unity, hence the decay into axion is negligible just at the 
end of thermal inflation. However, as the temperature decreases, the decay rate into 
the axion pair becomes dominant. In particular, even if T, ''^' becomes larger than the 
expansion rate just after the thermal inflation, the saxion coherent oscillation cannot be 
fully dissipated away because the dissipation rate becomes suppressed with the decrease of 
the cosmic temperature. (See Appendix.) Thus, the decay of the saxion may overproduce 
the axion. If the saxion still does not dominate the Universe when the perturbative decay 
into the axion becomes efficient {H ~ T^^2a), the axion overproduction is avoided. 

In order to estimate the present energy density of axion, we have performed numerical 
calculation. We have solved a set of equations 

p^ + i3H + rf ^^ + T^^2a)p^ = 0, (4.18) 

Prad + 4i7prad = rf ^V* + rf ^Va, (4.19) 

Pa + AH Pa = T^^2aP4> - rf ^Va, (4.20) 



3H'M^, = p^ + prad + Pa. (4.21) 

with initial conditions T = Tpx, Pa = and p^ = 'ni? fl, where ri ^^' denotes the axion 
dissipation rate, which is roughly sam( 
number of neutrino species AiVeff by 



dissipation rate, which is roughly same as T\ . Then we have calculated the extra effective 
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AiVeff = f ^^ M^ , (4.22) 

where g^:s is the effective number of massless degrees of freedom. The results are plotted 
in Figs. [21 In this plot we have taken fa/M = 100 (hence rris ~ 100m) and A = 1 (top) 
and A = 0.05 (bottom), and g^s = 228.75. The recent Planck result |46j excludes the 
parameter regions with AN^s ^1- It is seen that parameters with m > a few TeV and 
fa ^ 10^ GeV is viable, (/a ^ 4 x 10^ GeV is required from astr ophy sical arguments |22J.) 



Therefore, heavy SUSY scenario with m ^ 0(1) TeV is favored l^^i We have checked that 
the figure looks quite similar for fa/M = 1. To see more detail, we have plotted time 
evolution of various quantities as a function of Hpt/H in Fig. |3l Parameters are chosen as 
A = 1, m = 5 TeV, fa = 10^ GeV, M = 10'^ fa with which we obtain AN^s ~ 1- Since the 
dissipation rate T\ '^ is larger than the Hubble rate at early time, a significant amount of 



'^^'* Thermal inflation occurs for most of the parameter regions in the top panel of Fig. [5] In the bottom 
panel, thermal inflation does not take place in the about upper half of the parameter region (see Eq. (|4.3p ). 
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the saxion energy goes to the radiation. The produced axions are also thermahzed since 
ri ^^' is also large. Although gradually T\^ / H decreases, the perturbative decay T^^2a 
becomes efficient before the saxion dominates the Universe and hence the overproduction 
of axion can be avoided. 

Note that in a case with large dissipation rate T \ ^ if, axinos and saxions will 
be thermally populated. This does not cause cosmological problems as long as the axino 
is heavy enough to decay well before BBN and the produced LSPs have relatively large 
annihilation cross section |171 HHl IH]. All these necessary properties are consistent with 
heavy SUSY scenario, such as the pure gravity- mediation [50], which leads to 125 GeV 
Higgs pLj. 

5 Conclusions and Discussion 

In this paper, we have investigated the dynamics of a scalar field in thermal environment 
with a symmetry breaking potential. In particular, we have performed a detailed study of 
the trapping of scalar field at an enhanced symmetry point. 

Throughout this paper, we have focused on the case where the scalar field interacts with 
fermion and boson which are charged under SM gauge group and the initial amplitude of 
scalar field is far from the minimum of the potential. We considered the trapping dynamics 
of scalar field with the interaction of SM plasma taken into account. Although we have 
studied the SUSY PQ model in detail to make our discussion concrete, we emphasize that 
our results are rather general and that they are applicable to other symmetry breaking 
dynamics than the saxion. Since the trapping occurs due to the particle production, there 
are two possible sources: the production from the background thermal environment or the 
non-perturbative production from saxion itself. It is also notable that the higher order 
correction to thermal potential becomes important, so-called thermal logarithmic poten- 
tial. We have taken into account all these effects carefully and found that the dynamics 
of scalar field complicatedly depends on the interactions of produced particles to the SM 
particles. However, even if the saxion has a large initial amplitude, it is found that the 
saxion is likely trapped at its origin once and often leads to the thermal inflation for the 
most parameters of our interest. 

We have also studied the dynamics of saxion after the phase transition. It is noticeable 
that, even if the saxion once dominates the Universe and the thermal inflation occurs, 
the saxion can successfully dissipate its energy by the interaction with the thermal bath. 
In such a case, we can avoid the overproduction of the axion by the decay of the saxion 
field even if the dominant decay mode of the free saxion is into axion pair. To verify 
this statement, we have performed numerical calculation, and it is shown that the axion 
overproduction does not occur in the parameter region that are consistent with the high- 
scale SUSY scenario. 
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Figure 2: Contours of A^eff on {fa,m) plane. We have taken fa/M 
100m) and A = 1 (top) and A = 0.05 (bottom). 



100 (hence rris 
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Figure 3: Time evolution of various quantities as a function of Hp^/H 



normalized by the initial saxion energy density P0(ini), and r0_j.2a, T 
by H. Parameters are chosen as A = 1, rra = 5TeV, fa 



(dis) 

10^ GeV, M = 10-2/a. 



0) Prad) Pa 

Ta ^^ normalized 
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Appendix 

A Time Evolution of p^ 



in the presence of time-dependent dissipation effect 
0. (A.l) 



Let us consider the evolution of the 

r(t): 

Generally, r(t) is parametrized as r(t) = rj(tj/t)" where tj is an initial time. In the case 
of standard perturbative decay, n = 0. The case of n = 1 includes the Hubble friction: 
r(t) = Titi/t with Titi = 1/2(2/3) corresponding to the RD (MD) Universe. Thermal 
dissipation effect may correspond to n > 1. This equation can be easily integrated to 
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obtain 






t 

exp 



r,ti 




n-l 



1 



for n = 1, 
for n ^ 1. 



(A.2) 



From this expression it is easily checked that p^ exponentially decays with time for n < 1. 
For n > 1, the energy density exponentially decreases at t > (a few) x U if Fjtj ^ 1. 
Therefore, if Fjtj ^1,0 loses most of its energy within a few Hubble time. Eventually p<^ 
approaches to the asymptotic value 



p<t>{t) -> p^(ti)exp 
at t > tf where tf is defined by tfT(tf) = 1. 

References 



J- i^i 



n — 1 



(A.3) 



[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv:1207.7214| 
[hep-ex]]; S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) 
larXiv: 1207. 72351 [hep-ex]]. 

[2] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977); Phys. Rev. D 16, 
1791 (1977). 

[3] D. A. Kirzhnits and A. D. Linde, Annals Phys. 101, 195 (1976). 

[4] L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 76, 1011 (1996) 
|hep-th/9510119j; I. L Tkachev, Phys. Lett. B 376, 35 (1996) |hep-th795i0146j; 
S. Khlebnikov, L. Kofman, A. D. Linde and I. Tkachev, Phys. Rev. Lett. 81, 2012 
(1998) [hep-ph/980 4425]; G . N. Felder, L. Kofman, A. D. Linde and L Tkachev, JHEP 
0008, 010 (2000) | hep-ph/000402 4|. 

[5] L. Kofman, A. D. Linde, X. Liu, A. Maloney, L. McAllister and E. Silverstein, JHEP 
0405, 030 (2004) | hep-th/040300 1|. 

[6] K. Yamamoto, Phys. Lett. B 168, 341 (1986); G. Lazarides, C. Panagiotakopoulos 
and Q. Shafi, Phys. Rev. Lett. 56, 557 (1986). 



[7] D. H. Lyth and E. D. Stewart, Phys. Rev. Lett. 75, 201 (1995) |hep-ph/95024l7l . 
[8] D. H. Lyth and E. D. Stewart, Phys. Rev. D 53, 1784 (1996) |hep-ph/9510204 |. 



[9] A. Anisimov and M. Dine, Nucl. Phys. B 619, 729 (2001) |hep-ph/0008058l 



28 



[10] A. Berera, Phys. Rev. Lett. 75, 3218 (1995) |astro-ph/950"9049l ; A. Berera, I. G. Moss 
and R. O. Ramos, Rept. Prog. Phys. 72, 026901 (2009) [arXiv:'08083855! [hep-ph]]; 
M. Bastero-Gil and A. Berera, Int. J. Mod. Phys. A 24, 2207 (2009) |arXiv:0902.052T] 
[hep-ph]]. 



[11 
[12 

[13 

[14 
[15 

[16 

[17 

[is: 

[19 

[20 

[21 

[22 
[23 
[24 
[25 
[26 
[27 



J. 'i. Yokoyama, Phys. Lett. B 635, 66 (2006) |hep-ph/05 100 91]. 

M. Drewes, l arXiv:1012.5380l [hep-th]. 

M. Bastero-Gil, A. Berera and R. O. Ramos, JCAP 1109, 033 (2011) [arXiv:1008.1929i 
[hep-ph]]; M. Bastero-Gil, A. Berera, R. O. Ramos and J. G. Rosa, JCAP 1301, 016 
(2013) jarXiv:1207.0445] [hep-ph]]. 



T. Moroi and M. Takimoto, Phys. Lett. B 718, 105 (2012) [a rXiv: 1207.4858 ' [hep-ph]]. 

K. Mukaida and K. Nakayama, JCAP 1301, 017 (2013) |arXiv:1208.3"399l [hep-ph]]; 
JCAP 1303, 002 (2013) [a rXiv:1212.4985i [hep-ph]]. 

L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 73, 3195 (1994) 
|hep-th/ 9405187|. 

L. Kofman, A. D . Linde and A. A. Starobinsky, Phys. Rev. D 56, 3258 (1997) 
|hep-ph/9704452]. 



K. Enqvist, D. G. Figueroa and R. N. Lerner, JCAP 1301, 040 (2013) 
| arXiv:1211.5028l [astro-ph.CO]]. 

L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

G. N. Felder, L. Kofman and A. D. Linde, Phys. Rev. D 59, 123523 (1999) 
|hep-ph/98 12289] . 

For reviews, see J. E. Kim, Phys. Rept. 150, 1 (1987); J. E. Kim and G. Carosi, Rev. 
Mod. Phys. 82, 557 (2010) jarXiv:0807.3T25l [hep-ph]]. 

M. Kawasaki and K. Nakayama. larXiv:1301. 11231 [hep-ph]. 

T. Kugo, I. Ojima and T. Yanagida, Phys. Lett. B 135, 402 (1984). 

K. Tamvakis and D. Wyler, Phys. Lett. B 112, 451 (1982). 

J. E. Kim, A. Masiero and D. V. Nanopoulos, Phys. Lett. B 139, 346 (1984). 

K. Rajagopal, M. S. Turner and F. Wilczek, Nucl. Phys. B 358, 447 (1991). 

T. Asaka and M. Yamaguchi, Ph ys. Lett. B 437, 51 (1998) |hep-ph/9805449| ; Phys. 
Rev. D 59, 125003 (1999) |hep-ph/9811451| . 



29 



[28] N. Abe, T. Moroi and M. Yamaguchi, JHEP 0201, 010 (2002) [ hep-ph/OlllTSSl . 

[29] S. Nakamura, K. -i. Okumura and M. Yamaguchi, Phys. Rev. D 77, 115027 (2008) 
[ arXiv:0803.3725l [hep-ph] ] . 

[30] K. S. Jeong and M. Yamaguchi, JHEP 1107, 124 (2011) |arXiv:1102.3"30n [hep-ph]]. 

[31] K. Choi, E. J. Chun and J. E. Kim, Phys. Lett. B 403, 209 (1997) |hep-ph/ 9608"222|. 

[32] E. J. Chun, D . ComeUi and D. H. Lyth, Phys. Rev. D 62, 095013 (2000) 
|i hep-ph/0008133j . 

[33] S. Kim, W. -I. Park and E. D. Stewart, JHEP 0901, 015 (2009) |arXiv:0807.3607l 
[hep-ph]]. 

[34] K. Choi, K. S. Jeong, W. -I. Park and C. S. Shin, JCAP 0911, 018 (2009) 
[arX iv:0908.215 4 [hep-ph]]. 

[35] W. -I. Park, JHEP 1007, 085 (2010) [arXiv:1004.2326l [hep-ph]]. 

[36] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B 104, 199 (1981); A. R. Zhit- 
nitsky, Sov. J. NucL Phys. 31, 260 (1980) [Yad. Fiz. 31, 497 (1980)]. 

[37] M. Kawasaki, N. Kitajima and K. Nakayama, Phys. Rev. D 82, 123531 (2010) 
|arXiv:1008.5"0T3l [hep-ph]]: Phys. Rev. D 83, 123521 (2011) |arXiv:1104.1262i [Hep- 

ph]]. 

[38] N. Arkani-Hamed, G. F. Giudice, M. A. Luty and R. Rattazzi, Phys. Rev. D 58, 
115005 (1998) |hep-ph /9803290]. 

[39] K. Choi, E. J. Chun, H. D. Kim, W. L Park and C. S. Shin, Phys. Rev. D 83, 123503 
(2011) jarXiv:1102.2900l [hep-ph]]. 

[40] K. Nakayama and N. Yokozaki, JHEP 1211, 158 (2012) |arXiv:1204.5"420l [hep-ph]]. 

[41] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shifman, A. L Vainshtein, V. L Za- 
kharov, NucL Phys. B 166, 493 (1980). 

[42] S. Kasuya, M. Kawasaki and T. Yanagida, Phys. Lett. B 409, 94 (1997) 
|hep-ph/9608405]; Phys. Lett. B 415, 117 (1997) |hep-ph/9709202] . 

[43] S. Kasuya and M. Kawasa ki, Phys. Rev. D 56, 7597 (1997) [hep-ph/9703354] ; Phys. 
Rev. D 61, 083510 (2000) |hep-ph/9903324| . 

[44] L Tkachev, S. Khlebnikov, L. Kofman and A. D. Linde, Phys. Lett. B 440, 262 (1998) 
|hep-ph /9805209]. 

[45] M. Laine, Prog. Theor. Phys. SuppL 186 (2010) 404 |arXiv:1007.2"590l [hep-ph]]. 

30 



[46] P. A. R. Ade et al. [Planck Collaboration] , |ar2£ivil303,5082 [astro-ph.CO]. 

[47] H. Baer, A. Lessa, S. Rajagopalan and W. Sreethawong, JCAP 1106, 031 (2011) 
[arX iv:1103.54T3l [hep-pli]]. 

[48] H. Baer, A. Lessa and W. Sreethawong, JCAP 1201, 036 (2012) |arXiv:1110.2"49T] 
[hep-ph]]. 

[49] T. Moroi, M. Nagai and M. Takimoto. ErXiv:1303.0948l [hep-ph]. 

[50] M. Ibe and T. T. Yanagida, Phys. Lett. B 709, 374 (2012) |arXiv:1112.2"462] [hep- 
ph]]; M. Ibe, S. Matsumoto and T. T. Yanagida, Phys. Rev. D 85, 095011 (2012) 
|arXiv: 1202.22531 [hep-ph] ] . 



31 



